AN AUTOMATIC MESH GENERATOR 
USING TWO AND THREE-DIMENSIONAL 
ISOPARAMETRIC FINITE ELEMENTS 



James R. Adamek 



LIBIA'S— - 

NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIF. 93940 







L 



Monterey, California 



/ 




TH 



'It*— 





AN AUTOMATIC MESH GENERATOR 
USING TWO AND THREE-DIMENSIONAL 
ISOPARAMETRIC FINITE ELEMENTS 

by 



James R. Adamek 



Thesis Advisor: G. Cantin 



June 1973 




i . 



Approved ] {oh. pub-tic AcZcoac; diA&ubutlon uiilijnitcd. 



An Automatic Mesh Generator Using Two 
and Three-dimensional Isoparametric Finite Elements 



by 



James R. ^darnek 

Lieutenant, United States Navy 
B.S., University of Oklahoma, 1965 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1973 



LIBRARY 

na- t, :ostgh*buate schood 

1WIU .. , CALIF. 93940 



ABSTRACT 

The objective of the project described in this report 
was to develop computer systems which would generate the 
element connectivity, and nodal point co-ordinates for two 
and three-dimensional finite element programs using isopara- 
metric finite elements. The computer systems and sample 
problems are discussed. 
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INTRODUCTION 



I . 

The increased popularity of the Finite Element analyti- 
cal method due to the availability of high speed, large 
memory computers has led to the solution of many heretofore 
unsolvable problems in structural analysis. Although users 
of this method welcome the greater freedom given them to 
pursue more fundamental problems, they are presented the 
tedious task of generating and checking the mesh input data 
of their particular computer system. This data includes the 
arrangement of nodes within each finite element (connecti- 
vity) and the co-ordinates of all the nodes. 

Automatic mesh generation is an attempt at simplifying 
input data for the finite element programs. If the user 
can communicate with the computer in a better, faster, and 
more efficient manner, money and labor will be substantially 
saved, and the possibility of human error will be greatly 
reduced . 

Many automatic mesh generation schemes exist, however 
these systems are normally proprietary, very specialized, 
and not generally adaptable to finite elements having side 
nodes. The present system is an attempt to remedy that 
situation . 

References 3 and 6 prompted this author to devise a code 
for two and three-dimensional automatic mesh generation 
schemes, utilizing the technique of mapping from a local 
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non-dimensional system of reference into the global cartesian 
system of reference. 

Reference 4 is a previous work on a three-dimensional 
automatic mesh generation scheme using isoparametric elements, 
however the coding incorporated an interpolation scheme which 
led to ill-conditioned elements in some instances. 

To illustrate the effect of reduction in the amount of 
data as a result of automatic mesh generation, the following 
example is introduced: Figure 1(a) is a rather complex 

geometrical shape in the x-y plane, which can be sub-divided 
into sections (super-elements) for which one of the isopara- 
metric elements is a very close geometric approximation. 

These super-elements may be further discretized into smaller 
elements. The connection of the super-elements and elements 
' is more easily illustrated in a local non-dimensional system 
of co-ordinates, as shown in Figure 1(b), where all the super- 
elements are squares. With the appropriate boundary descrip- 
tion of the super elements as shown in Figure 1(a), and a 
super-element connectivity diagram as shown in Figure 1(b), 
the fine grid of figure 1(b) can be mapped back into the 
original object to produce the desired mesh of Figure 2. 

Four automatic mesh generation schemes will be presented 
in this thesis. One is for two-dimensional linear, quadratic 
and cubic elements. The others for three-dimensional linear, 
quadratic and cubic elements respectively. 
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Original object sub-divided into super-elements 

Ka) 




Connectivity of super-elements in local system 

Kb) 



Figure 1. SEL Diagrams. 




Figure 2. Final Mesh. 
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II. ISOPARAMETRIC FINITE ELEMENTS 



The isoparametric concept will be discussed in this 
chapter in sufficient detail to establish only the procedures 
for constructing isoparametric elements and the validity of 
relationships used to obtain a mapping from the local non- 
dimensional reference into the global cartesian system of 
reference. The readers are referred to references 1 and 2 
for further details of the isoparametric concept. 

The serendipity family of isoparametric elements [1] was 
utilized in this paper due to the ease of transformation 
(mapping) from a local non-dimensional system of coordinates, 
Avhere all elements are squares or cubes, to a cartesian 
system with curvilinear boundaries. 

A. TWO-DIMENSIONAL ELEMENTS 

Consider the square isoparametric elements of Figure 3. 
Figure 3(a) is a linear element consisting of four corner nodes. 




linear element 
3(a) 





3(b) 



cubic element 
3 (c) 



Figure 3. Isoparametric Element Configurations 
in Local (£,q) Coordinates. 



Figure 3(b) is a quadratic element consisting of four corner 
nodes and four side nodes. Figure 3(c) is a cubic element 
consisting of four corner nodes and eight side nodes. 

A curvilinear element, such as that of Figure 4, may be 
produced by mapping the appropriate square of Figure 3 into 
the x-y plane using the relations 



x 



n 

Z N.x. 



i= 1 



i i 



n 

y = Z 
i=l 



N . y . 

i' i 



CD 



where (x^,y^) are the cartesian coordinates of node i. 
hh are shape functions of the serendipity family, listed in 
Section II.C.l., n is the number of nodes on the element 
(i.e., 4 for linear elements, 8 for quadratic elements, 12 
for cubic elements) . 

Note that the convention for node numbering within any 
element (see Figures 3 and 4) is to begin at (£»n) = (1,1) 
and to number the nodes proceeding counterclockwise. 




Figure 4. Curvilinear Element Mapped 
in Cartesian Coordinates. 
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B. THREE-DIMENSIONAL ELEMENTS 



Linear elements. Figure S, have eight corner nodes. 
Quadratic elements, Figure 6, have eight corner nodes and 
twelve side nodes. Cubic elements, Figure 7, have eight 
corner nodes and twenty-four side nodes. 

Similarly, as with two-dimensional elements, three- 
dimensional curvilinear elements may be produced (mapped) 
using the relations 



x 



n 

Z 

i=l 



N i*i 



n 

y = Z 
i = l 



Vi 



n 

z = Z 
i=l 



N i z i 



where (x^,y^,z^) are the cartesian coordinates of node i. 

As shown in Figures 5, 6 and 7, the convention for node 
numbering within any element is to begin at (£,n,C) = (1,1,1) 
and to number the nodes proceeding counterclockwise about 
the zeta axis. 

6 (- 1 » 1 > - 1 ) 




15 




Figure 6. Three-dimensional Quadratic Element. 
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Figure 7. Three-dimensional Cubic Element. 
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C. SHAPE FUNCTIONS 



Let £ 



o 



= «i 



= nn. 



’o 






1 . Two-Dimensional Elements 

a. Linear elements 
♦ 

See Figure 3(a) for identification of the nodes. 
N i = (1/4) U o + 1) (n Q +1) (3) 



where = ±1 and n- = ±1. 

'i 

b. Quadratic elements 

See Figure 3(b) for identification of the nodes. 
Corner nodes: 1,3, 5, 7 

N i = (1/4) (l+? 0 ) (l + n Q ) (C 0 + n 0 - 1) (4) 

where £ - = ±1 and n- = ±1. 

l l 

Side nodes: 2, 4, 6, 8 

Nj = (1/2) (1 - £ 2 ) (1 + n o ) , (5) 

for nodes 2 and 6, .where n ^ = ±1 

Ni = (1/2) (1 + C 0 ) (1 - r, 2 ) , (6) 

for nodes 4 and 8, where ^ = ±1. 

c. Cubic elements 

See Figure 3(c) for identification of nodes. 
Corner nodes: 1,4,7,10 

N t = (l/32)(l + ? o )(l+n o )[9(C 2+ n 2 ) - 10] (7) 

where ^ = ±1, rii = ±l. 
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Side nodes: 5,6,11,12 



N i = (9/32)(l+£ 0 )(l-n 2 )(l+9n 0 ), (8) 

where = ±1, n ^ = ±1/3. 

Side nodes: 2, 3, 8, 9 

= (9/32)(l+n o )(l-? 2 )(l + 9? 0 ), (9) 

where E,. = ±1/3, n • = ±1. 

1 * 1 

2 . Three-dimensional Elements 

a. Linear elements 

See Figure 5 for identification of nodes. 

Ni = (1/8) (1 + £ q ) (l+n Q ) (1 + C 0 ) , (10) 

where = ±1, r) ± = ±1, = ±1. 

b. Quadratic elements 

See Figure 6 for identification of nodes. 
Corner nodes: 1,3,5,7,13,15,17,19 

Ni '= (1/8) (in o ) Ci+n 0 ) dn 0 ) (5 0 + VV 2 > 

where = ±1, r\ i = ±1, = ±1. 

Side nodes: 2,6,14,18 

Ni = (1/4) (l-£ 2 ) (l + n 0 ) (1 + C 0 ) (12) 

where £ ± = 0, r\ i = ±1, ^ = ±1. 

Side nodes: 4,8,16,20 

Ni = (l/4)(l-n 2 )(m 0 )d + C 0 ) (13) 



( 11 ) 
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where ^ = ±1, n i = 0, c i = ±1 
Side nodes: 9,10,11,12 

N i = (1/4) (1-c 2 ) (l + c 0 ) (l + n 0 ) 
where £. = ±1, ri . = ±1, C- = 0 
c. Cubic elements 

See Figure 7 for identification of nodes. 
Corner nodes: 1,4,7,10,21,24,27,30 

N i = (1/64) (1H 0 ) (l + n o ) (1 + C Q ) [9 (5 2 + n 2 + C 

where = ±1, n- = ±1, = ±1 

Side nodes: 2,3,8,9,22,23,28,29 

N i = (9/64)(l-C 2 )(l+9C 0 )(l+n o )(m o ) 

where £ i = ±1/3, r\ i = ±1, ^ = ±1 
Side nodes: 5,6,11,12,25,26,31,32 

N i = (9/64)(l-n 2 )(l + 9n 0 )(l + C 0 )d + C 0 ) 

where = ±1, n i = ±1/3, = ±1 

Side nodes: 13,14,15,16,17,18,19,20 

N i = (9/64)(l-c 2 )(l + 9 ?o )(l + C 0 )(l + ri 0 ) 
where £ i = ±1, = ±1, K i = ±1/3. 



( 14 ) 



2 ) - 19] 
(15) 



(16) 



(17) 



( 18 ) 



III. BASIS OF AUTOMATIC MESH GENERATION SCHEME 



As mentioned in the introduction, reference 3 prompted 
the utilization of the following scheme in this thesis 
research . 

The essence of the scheme is the use of simple mapping 
functions to transform simple orthogonal connectivity dia- 
grams into real geometrical shapes. Considering the partic- 
ular two-dimensional case of a parabolic quadrilateral of 
Figure 4, in which the x and y co-ordinates of eight nodes 
are known, we can write equation Cl) from Chapter II, where 
N^ is a shape function associated with each node and defined 
in terms of a local non-dimensional coordinate system £ and n 
which has values ranging from 1 to -1 on opposite sides. 

Typical shape functions are given in II. C. , and for this 
example equations (4), (5), and (6) apply. If the co-ordinates 

of the nodal points are known, then cartesian co-ordinates of 
any specified point (£,n) can be found by equation (1). 

If the whole region in which the mesh is to be generated 
could be described adequately by a quadrilateral of the shape 
in Figure 4, a mesh of any refinement could be automatically 
generated inside it by specifying the coordinates of the 
eight boundary nodal points and the number of sub-divisions 
in the £ and n directions. The above scheme is shown pic- 
torially in Figures 8(a) and 8(b). 
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► x 




Quadratic Super-Element 
local coordinate 
system 

8(a) 



Quadratic Super-Element 
cartesian coordinate 
system 

8(b) 




8(c) 



Super-Element Type 



A 

B 

C 

D 

E 



quadratic 

linear 

linear 

cubic 

quadratic 



A corresponds to corner nodes 
© corresponds to side nodes 

[) corresponds to side nodes for a single Super-Element 
whose adjacent Super -Element is of a different type 



Figure 8 



17 



■mMI 













If the whole region could not be adequately described by 
Figure 4, then we could subdivide the region into more than 
one quadrilateral and apply the same principles to each as 
specified above. These region subdivisions will be referred 
to as Super-Elements hereafter (see Figure 8) . 

The local Super-Element of Figure 8(a) can be quite 
drastically distorted, even up to a point of making two sides 
lie along the same line in the cartesian system such as in a 
triangle or semi-circle. Care must however be taken not to 
make any corner angles greater than 180° as a non-uniqueness 
may result [1]. In addition, side nodes should be equally 
spaced between corner nodes to insure equal spacing of sub- 
divisions in the £ and n directions and thus the x and y 
directions in the mapped cartesian system. 

The criteria above also apply to the linear and cubic 
isoparametric Super-Elements for two-dimensional meshes, and 
the linear ,' quadratic and cubic isoparametric Super-Elements 
for three-dimensional meshes. 

The decision on whether to use linear, quadratic or 
cubic Super-Elements will depend entirely upon the region 
boundary to be mapped. Linear Super-Elements will map 
quadrilaterals with straight sides; quadratic Super-Elements 
will map quadrilaterals with straight or parabolic sides; 
cubic Super -El ements will map quadrilaterals with straight, 
parabolic or cubic sides. Figure 8(c) is an illustration of 
Super-Element selection based upon region (structure) boun- 
dary conditions. 
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The final generated mesh for a given structure may con- 
tain nodes on the boundary which do not exactly coincide 
with the actual structure boundary, or the elements within 
may not be represented sufficiently accurately. In this 
case it may be necessary to have some of the coordinates 
adjusted once the mesh has been generated. This adjustment 
will be quite easily accomplished with the aid of computer 
plots, which can be obtained as output to the mesh generat- 
ing programs discussed in Chapter IV. 
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IV. DISCUSSION OF COMPUTER PROGRAMS 



The computer programs presented in this thesis were 
written to support the computer systems 'PLISOP' and 'TRISOP', 
both of which were coded by Professor G. Cantin, and his 
students of the Naval Postgraduate School, Monterey, Cali- 
fornia. TRISOP performs a structural analysis of three- 
\ 

dimensiona.l problems using isoparametric finite elements, 
while 'PLISOP' performs the same analysis of two-dimensional 
problems. The input data to each system is quite extensive 
with regard to element connectivity and nodal point coordi- 
nates . 

The automatic mesh generating systems discussed in this 
chapter compute the element connectivity and nodal point 
coordinates required for the above two computer systems. 

Four mesh generating programs will be discussed. The first 
generates data for 'PLISOP' with linear, quadratic and cubic 
elements, and is called 'PLIMEG 1 . The other three programs 
generate data for 'TRISOP' with linear, quadratic, and cubic 
elements respectively. They are called 'TRIMEG-1', 'TRIMEG-2' 
and 'TRIMEG-3'. 

A. THE GENERAL PROGRAM 

While a simple mesh generator using only one super- 
element may be sufficient for some purposes, it is limited to 
simple geometries with but a single uniform property. For 
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complete generality, as suggested at the end of Chapter III, 
the scheme was extended to many super-elements, either linear, 
quadratic, cubic or any combination of the three, and repre- 
sented in a 'checkerboard' pattern, where each of the super- 
elements may be sub-divided and each may define a material 
with a single property. The above scheme is presented in 
Sections B and C of this chapter for two and three-dimensional 
meshes respectively. 

B. TWO-DIMENSIONAL MESH GENERATION (PLIMEG) 

1. Definitions (see Figure 9) 

a. SEL: an acronym, used to refer to a super- 

element as defined previously on page 7. 

b. Master Grid: A 10 x 10 'checker-board' array 

of SELs which is utilized as a working 
board for SEL numbering and orientation. 

The master grid SELs are numbered sequen- 
tially along grid columns. 

c. Rows: subdivisions of a SEL in q direction 

d. .Columns: subdivisions of a SEL in £ direction 

e. Elements: those portions of a SEL created by 

rows and columns 

2 . Super-element Selection 

The number of SELs required in a given mesh is deter- 
mined from boundary considerations and the types of materials 
from which the object is composed. Enough SELs must be 
incorporated to adequately define the boundary of the object, 
and at least one SEL must be utilized for each different type 
of material. Figure 10 shows an example for a particular 
object. The reader will note in Figure 10(a) that the 
common boundary of SEL 1 and SEL 2 is based on material type, 
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SEL #12 (8 nodal point) 
9(b) 



^indicates corner nodes; 0 indicates side nodes 
All nodes indicated comprise SEL boundary nodes- 



Figure 9. PLIMEG Definitions. 



while the remaining boundaries are determined from geometri- 
cal considerations only. 

Once the quantity and type of SELs have been deter- 
mined, the proper SEL orientation in the local nondimens ional 
system of reference must be decided upon, as shown in Figure 
10(b). The local system is now superimposed upon the master 
grid of Figure 9 to obtain the SEL numbering scheme of 
Figure 10(b). Note: this superimposition may be anywhere 

within the boundary of the master-grid. 

3 . Mesh Criteria 

The following rules must be adhered to at all times: 
a. All elements within a given mesh must be of the 
same type (i.e., linear, or quadratic or cubic), and only 
the type of SEL may vary. 
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A corresponds to corner nodes 

G corresponds to side nodes 

ID corresponds to side nodes for 
a single SEL whose adjacent 
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Local SEL Diagram 
10(b) 




Figure 10. PLIMEG Example. 
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b. The number of rows in all SELs of any master grid 
row must be equal (i.e., SELs 2 , 12 , 22 . . . 92) . 

c. The number of columns in all SELs of any given 
master grid column must be equal (i.e., SELs 11 , 12 , 13 . . . 20) . 

d. SELs may be connected to adjacent SELs in the 

following ways (see Figure 11) : totally connected, connected 

at one corner only, totally disconnected, and connected at 
the corner nodes only but not along common edge. Figure 11 
shows only horizontally connected SELs, but the same connec- 
tion schemes apply to vertically connected SELs as well. 

4 . Bandwidth, Node and Element Numbering 

Element numbering is always along master grid 
columns, and with the exception of discontinuities between 
adjacent SELs, such as in Figures 1 1 Ct> ) and 11(d), node 




Totally 

Connected 



Partially Totally 

Connected Disconnected 



11(a) 



11(b) 



11(c) 



V 

V 




Connected at 
Corner Nodes Only 
11(d) 



Figure 11. SEL Connections. 
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numbering is also along master grid columns. Therefore, to 
obtain a node numbering scheme which results in the best 
computational efficiency, since 'PLISOP' and 'TRISOP' are 
designed to be used with a banded stiffness matrix solution, 
it is prudent that the number of rows in the master grid row 
containing the most columns be less than or equal to the 
total number of columns in tfrat particular master grid row. 
(The reader is referred to Section C.3 of this chapter for 
additional information regarding this topic.) For example, 
in Figure 10 the master grid row containing the most columns 
is grid row two, in which there are sixteen columns, while 
the number of rows in this grid row is five, thereby satis- 
fying the above criteria for bandwidth. Another example 
concerning bandwidth is as follows: suppose the object to 

be considered is in the shape of an "L" and that all SELs 
contain an equal number of rows and columns (i.e., 3). Two 
different SEL combinations could be utilized to describe 
the mesh as shown in Figures 12(a) and 12(b). However it is 
evident from Figures 12(c) and 12(d) that bandwidth would be 
greater for the shape of Figure 12(b), since there are six 
rows in master grid column one for Figure 12(d) vice three for 
Figure 12(c) . 

There will be many instances when the geometries of 
Figure 12 will present itself, and when this condition arises, 
the user must determine whether element shape within the mesh 
or a mesh with minimum bandwidth is the more crucial. 
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Local SEL Diagram 
12(c) 
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Cartesian SEL Diagram 

12(b) 
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i 



Local SEL Diagram 
12(d) 



Figure 12. The "L" Shape. 



Whenever discontinuities such as those of Figures 
11(b) and 11(d) exist in a mesh, a good guideline to follow 
is to have them exist above or below a SEL, if possible, 
rather than to the left or right, since the latter type of 
SEL connection occasionally results in a somewhat non- 
efficient node numbering scheme due to the complicated 
coding incorporated in Subroutine CONN of PLIMEG. 

5 . PLIMEG Output 

The major portion of the output with PLIMEG is 
element connectivity and nodal coordinates, however addi- 
tional output has been incorporated to aid in the use of 
PLISOP. To begin with, all data input to PLIMEG will be 
reproduced as output as a check to the user, in event a 
mispunched data card is experienced. The total number of 
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elements and joints in the mesh together with the half 
bandwidth will also be printed. In addition, an option to 
punch on data cards the element connectivity and coordinates 
of the joints is incorporated. Lastly, a plot of the mesh 
may be obtained on the printer, whereby Subroutine GRID calls 
Subroutine UTPLOT in the IBM 360 source library. However, 
since plots on the printer are somewhat distorted and only 
consist of points, a plot can also be obtained on the off- 
line plotter (Calcomp) whereby Subroutine GRID calls Sub- 
routine DRAW in the source library. The resulting plot is 
as those shown in Figures 2 and 10(c). The printer plot 
should be used in determining proper mesh discretization and 
then Calcomp for the final plot of the mesh. Both plots can 
not be obtained on the same computer run. 

Appendix A discusses input data preparation for 
PLIMEG while Appendix B discusses sample problems and shows 
the input/output formats, the node and element numbering 
schemes, and the resulting plots. These two Appendices may 
be removed and utilized as an instruction manual for PLIMEG 
once the contents of this and the previous chapters are 
comprehended by the user. 

C. THREE-DIMENSIONAL MESH GENERATION (TRIMEG) 

TRIMEG actually consists of three separate programs: one 

for each type of element (i.e., linear, quadratic, or cubic). 
They are called TRIMEG-1, TRIMEG-2, and TRIMEG-3 respectively. 
However, since the computer scheme and the data input incor- 
porated in these programs is identical, they will be discussed 
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as one (i.e., TRIMEG). TRIMEG also incorporates many of 
the ideas of PLIMEG, therefore only the dissimilarities 
between the two schemes will be presented in this section. 

1. Definitions (see Figure 13) 

a. Master Grid: a 5x5x5 three-dimensional "checker- 

board" array of SELs. The master grid SELs are numbered 
sequentially along master grid columns in successive master 
grid slices . 

b. Slices: sub-divisions of SEL in i ; direction 

c. Element: those portions of a SEL created by 

rows, columns and slices 

2 . Mesh Criteria 

The following rules mus t be adhered to at all times. 

a. The number of rows in all SELs of any master grid 
row must be equal (i.e., 1 , 6 . . . 21 , 26 , 31 . . . 46 , 51 , 56 . . . 71 , 76 , 81 , 
.. .96,101,106, . . .121) . 

b. The number of columns in all SELs of any master 
grid column must be equal (i.e., 1-5,26-30,51-55,76-80, 
101-105). 

c. The number of slices in all SELs of any master 
grid slice must be equal (i.e., 26 , 27 , 28 . . . 50) . 

d. The faces of adjacent SELs may be connected in 
any of the ways shown in Figure 11 (applies to all six faces) , 
with the exception of Figure 11(d). If this situation exists 
in a mesh, then four SELs must be utilized as shown in 
Figure 14. 
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SEL #7 (20 nodal points) 
13(b) 

Figure 13. TRIMEG Definitions. 
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Figure 14. TRIMEG SEL Connection for 
Similar Connection of 
Figure 11 (a) . 

3 . Bandwidth, Node and Element Numbering 

Node and element numbering is identical to that of 
PLIMEG for each successive slice in the mesh. The following 
rules apply to obtain minimum bandwidth. The number of rows 
in the master grid row containing the most columns should be 
less than or equal to the total number of columns in that 
particular master grid row. In addition, the number of 
columns in this master grid row should be less than or equal 
to the total number of slices in this particular grid row. 
This scheme can be visualized by sketching a single SEL with 
one row, two columns and three slices and then numbering the 
nodes along columns in successive slices. The half-band 
width is 30 as shown in Figure 15(a). Using two rows, one 
column and three slices, produces a half bandwidth of 33, 
etc. The reader should sketch several possibilities for 
row, column and slice combination (Figure 15) to convince 
himself that orientation of the mesh determines bandwidth, 
and thus efficiency, when solving problems with PLISOP and 
TRISOP. As can be seen in Figure 15, numbering is identical 
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half-bandwidth = 30 
15(a) 



half-bandwidth = 33 
15(b) 
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3x2x1 

hal f-bandwidth = 54 

15(c) 



1x2x3 

half-bandwidth = 36 
15(d) 



Figure 15 
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to that of PLIMEG for a single SEL with one slice without 
the third dimension. For PLIMEG the half bandwidth for 
Figure 15(a) would be 8. 

4 . TRIMEG Output 

TRIMEG output is identical to PLIMEG except for the 
plotting package. With TRIMEG a plot can only be obtained 
on the Calcomp plotter. The plots are of three different 
projections on the planes produced by rotating the mesh 
about the z-axis, the latest x-axis and the latest y-axis 
respectively. The three angles of rotation (Euler angles) 
and the resulting axes are shown in Figure 16 and the scheme 
is contained in subroutine TRFR of TRIMEG. 




Figure 16. Axes Rotations for Projection Plotting. 
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Appendix C discusses input data preparation for 
TRIMEG while Appendix D discusses sample problems and shows 
the input/output formats, the node and element numbering 
schemes, and the resulting plots. These two Appendices may 
be removed and utilized as an instruction manual for TRIMEG 
once the contents of this and the previous chapters are com- 
prehended by the user. 



v 
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V. CONCLUSIONS 



The utilization of PLIMEG and TRIMEG automatic mesh 
generation schemes along with existing finite element pro- 
grams, now makes it much more practical to solve realistic, 
complicated problems, since the amount of input data has 
been reduced as much as 95 per cent. This decrease in the 
amount of input data cards required for the finite element 
programs not only reduces the chance for error, but makes 
problem solving more of a challenge and not just a headache. 
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VI. RECOMMENDATIONS 



PLOT OF STRUCTURE MESH 

In addition to the present plotting routine incorporated 
in TRIMEG, it would be desirable to obtain a plot of the 
generated mesh which did not contain hidden lines. Profes- 
sor Gilles Cantin of the Naval Postgraduate School, Monterey, 

\ 

California is presently working on such a solution. Adding 
the actual node numbers to the computer produced plots would 
also be of great value. 



35 



APPENDIX A 



INPUT DATA PREPARATION FOR f P LI MEG f 





St 


ft ft ft ft 














2 LO 


23 








ft 


ft St 






• 








OlO> 


03 








St 


ft ft 






o 








333 < 


3 


• 






ft 


ft ft 






2 






• 


302 2 


303 


>- 






ft 


ft 3ft 






CL 






3 


*-H — « 3 


— 33 


3 






St 


ft Oft 






X53 






2 


-£<tO< 


22 — 


•2 






Si 


ft U ft 






uo 






3 


3 • 


— 2 


XO 


2 




ft 


ft Xft 






30 






X 


LOO ox 


L0< 


LO 


— 




ft 


ft Oft 






3 






3 


Q333L0 


3 1 — 3 


3K 


• 




st 


ft LOft 






CL< 






3 


33333 


030 


22 


OLO 




ft 


ft ft 






oo 






3 


Cfvl<x CJ2 


303 


< 


CCL 




St 


ft 3ft 






N 








O X 


3 3 


XCL 


03 




St 


ft 3 ft 


• « 




03 






3 


!£3L03 


3 


3 O 


33 




ft 


ft <ft 


LO 




23 






3 


000 X 


03< 


— <X 


2 




St 


ft 3ft 


j£ 




< 






3 


23 03 


2 


20 


-O 




ft 


ft Oft 


3 




LO 








3 CO 3 


3-* 


O 


— 2 




■ft 


- ft <1 ft 


3 




30 






LO 


- 3 3 X 


-3C\J 


O 


in 




ft CL ft CL ft 


3 




22 




• 


3 


X3hNh 


V M <w » 


33 


—3 




ft O ft 3 ft 


O 




<< 




LO 


2 


22< — < 


22 


— LO 


02 




ft LO ft ft 


3 




20 




3 


►— < 


<3Z 3 3T 


< LO 


33 


33 




ft 


3ft h- ft 










2 


O 


33 — 0 


33 — 


LO — 


—2 




ft 3ft LO ft 


LO 




L03 




3 


3 


330 Q 


30 


— < 3 


3 




ft CX ft 13 ft 


<1 




- 2 — 




2 




3 3 23 


33 


1 — 


33 




St 


- ft LL ft 






23 m 




3 


3 


3 3 LJ 2 — 


2 CL (3 • 


<12 


<X LU 




ft 


ft ft 


lo 




3 




3 


< 


O^O <3 


U 1 CL LO — 


2 




St 


2 ft 3ft 


— 


• 


003 O 




3 


O 


OO3L0 


<0.3 




23 




ft Oft < ft 




o 


30 CL 






a 


O CO — 


O 233Q 


03 




ft LL ft >ft 


51 


CL 


< 




DC 


2 


3 L0> 3 


CL— 3 — 3 


23 




ft 


ft <ft 


< 


< 


• 30 




3 




3— -X< 


333033 


O 




ft 23ft 2 St 


CL 


O 


X2 




3 


3 


rvl XOLO 


M— — 323 


-LO 


m 


ft < ft St 


o 




30 3 


• 


3 


O 


X 2 


300 


LO 


o 


ft CL ft -ft 


o 


X 


3 3 




LO 




LOOQ3LO 


LO LO • 


3 


03 


LU 


ft Oft (MSt 


CL 


o 


32 LO 


in 




3 


— 2 23— • 


— — 1 — 3 


— 3 


20 


53 


ft Oft 3 St 


CL 


<1 


OOw 




3 


3 


OC2 


32 


4 


< 


— 


ft CL ft O' ft 




3 


3 


-r 


O 


3 


LO 3 3 


LO! — OO 


— 3 


02 


3 


ft CL ft 3 ft 


LO 




30 • 


w 




2 


hZ> *2 


— O 30 


3 


3 


CL 


ft 


•it -it 


►— < 


— * 


3 Q 




3 


3 


X 300 


X3333 


30 




ft 


Oit aCit 


X 


CO 


303 


3 


IU 


2*M 


3 3—3 


3 3 3 <1 3 3 


32 




ft 2 ft 3 ft 


3 


<x 


030 


< 


3 


—4 


h>X3 


200 


3 


LOO 



ft -*St CO St 



o 



30 



LU 



LL hQm LL O — 



2 

LU 

CO 

2 



LJ 



< 

O 



LU 

Q 



LU 

2 

LU 

3 
LU 

2 

UJ 

3 

3 



ft 3 ft 2 ft 


2 


— 


303 


2 


O 


1 2 


— • F— 2 |— 


•“-*2 22 3 12 


0 


LO 


ft < 


Si- 3 si- 


0 




320 


0 


2 


3 0 


0002 


— 300 


22 




ft 2 ft O -Si- 


3 


3 


3— 2 


2 






• * 3330 


3 OLO 






ft 3ft 3ft 




< 


— — 






CO 


-»I2 


• • *20X3 


Q 


• 


ft 2ft Oft 


3 


2 


33 




• • 




X 02 30 


— 3H3L0O 


32 


O 


ft 3 


ft ft . 


O 


2 


02 


0 




•• *4* 


020 03 


3022 3 


<3 


2 


ft Oft -ft 


3 


O 


2 O 


2 


3 


. — . 


200 lo 


0< 3 — 3 3 


30 




ft 


it it 


2 


3 


33 — 


3 


3 


3 H 


03 30 


32 <0 


OLO 


3 


ft X ft 3 Si- 


— 




— L03 


3 


LO 


33 


3 23 


3 3333 


3 <£ 


3 


ft LO 


it Sit 




•* 


203 


2 


2 


22 


230— 3 


23X03 3 




LO 


it Ul it < it 


> 


3 


C02O 


3 




— 2 


>-*32 3: CO 


w0302< 


•* 




•it Sit Q-it 


2 


3 




2 










2 




ft 


ft <ft 


< 


3 




< 










O 




ft 


ft • st 


LO 


— 




2 










3 




ft 


Si- CL ft 


LO 


3 


CO 


< 


in 


O 


in 


O 


Q 


in 


ft 


Si- ‘Si- 


3 




M- 


3 






3 


c\J 






ft 


ft oft 


0 


*» 














3 




ft 


ft Si- 


3 


CM 


O 


X 


0 


O 


0 


O 


2 


0 


ft 


ft >-st 


2 




3 


LO 


3 


3 


3 


3 


3 


3 


Si- 


ft CD Si- 




O 




3 










2 




ft 


ft Si- 


3 


2 


—4 


2 


«— < 


O 


*— 4 


O 


3 


—4 


Si- 


Si- Oft 


X 


< 












^4 


3 




ft 


ft 3 St 


3 




• 


•* 


• 








3 


# 


ft 


St 0 Si- 




^■4 


3 


CO 


3 










3 


ft 


ft ost 






0 




0 








2 


0 


ft 


St ost 




0 


0 


Q 


0 








3 


0 


ft 


ft St 




2 




2 










3 




ft 


ft ft 




<i 




C 










O 




ft 


ft ft 




0 




O 










LO 





LU 

23 

LU 

3 

LU 

2 

LU 

CL 

3 

LO 



LO 

12 

o 

a: 

LL 

a 

CL 

LU 

CD 

2 

3 

2 



-2 

o 

CL 



* * * 
St St ft ft ft 



36 



(COL): NUMBER OF COLUMNS IN SUPER ELEMENT 
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APPENDIX B 



PLIMEG EXAMPLES, INPUT AND OUTPUT 

This appendix discusses two example problems. The first 
will show SEL selection and connection procedures dealing 
with an eight SEL object with discontinuities along some 
edges, and also the node and element numbering schemes incor- 
porated in PLIMEG. The second will be a more straightforward 
example illustrating the detail of input data and output 
formats . 

A. EXAMPLE ONE (Padeye; see Figures 17 and 18) 

The common boundary of SEL 21 and SEL 22 is based on 
material type, while the remaining SEL boundaries are deter- 
mined from geometrical considerations only. As noted in 
Figure 17(a), the common boundary of SEL 2 and SEL 42 must 
be connected manually once the mesh has been generated. 

This is accomplished as follows (see Figure 17(c)). The 
final mesh of Figure 18 will not contain nodes 77-80, there- 
fore remove these nodes from the output punched deck of 
coordinates and change the element connectivity cards for 
elements 55,56 and 57. In other words, node 77 becomes 
node 1, node 78 becomes node 2, node 79 becomes node 3 and 
node 80 becomes node 4. The number of nodal points (joints) 
will now be 76 vice 80, while the number of elements will 
remain the same (i.e., 56). Half bandwidth will be increased 
from 20 to 148, but may be held to 50 if SELs 1,2,3,4,5,12,13 
and 14, for example, are used to describe the mesh. 
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Material 2 
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Figure 18. Final Mesh for Example One. 

The remaining variables for input to PLIMEG will be shown 
in the more straightforward example which follows. 

B. EXAMPLE TWO (See Figures 19, 20 and 21) 

Figure 19 shows a simplified mesh generation scheme con- 
sisting of two SELs, six elements (linear) and 12 nodal 
points (joints). Although SELs connected vertically could 
have' been utilized, this would have established a mesh with 
three rows and two columns resulting in a half-bandwidth 
greater than that of the existing mesh. 

Although the mesh is quite coarse, it adequately describes 
the data input for the given boundary situation. The SELs 
could have been subdivided to produce a much finer mesh, 
however the amount of input data would remain unchanged. 

Figure 20 shows in detail the input and input format 
required for PLIMEG as described in Appendix A, while Figures 
21 and 22 present the output and output format. The last 
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two columns of the connectivity matrix (Figure 22) are KIND 
and TYPE. These quantities represent the kind of element 
and the type of material for that element, respectively. 
KIND equal to one implies linear elements, two implies 
quadratic, and three implies cubic. This information is 
required input for PLISOP. 
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Figure 20. PLIMEG Input Data (example two). 
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Figure 21. PLIMEG Output (input data check). 
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Figure 22. PLIMEG Output. 
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APPENDIX D 



TRIMEG EXAMPLES, INPUT AND OUTPUT 

Since the coricept for TRIMEG is very similar to that of 
PLIMEG, only the figures illustrating two examples will be 
shown in this appendix. The user should be completely 
familiar with the contents of Appendix B (PLIMEG Examples, 
Input and Output) prior to utilizing this appendix and 
Appendix C (TRIMEG Input Data Preparation) . 

A. EXAMPLE ONE (Crack Propagation Problem) 

Figure 23 shows the local SEL diagram indicating SEL 
numbers and SEL corner node numbering. The cartesian SEL 
diagram is not shown since it would be identical to that of 
the local diagram since the object is a cube. Figure 24 
shows the final mesh, where each SEL contains 2 rows, 2 
columns and 2 slices. 

Example two will completely illustrate all input, output 
and respective diagrams. 

B. EXAMPLE TWO (Folded Plate) 

Figures 25, 26 and 27 show a simplified mesh generation 
scheme consisting of three SELs, 27 elements (linear) and 
eighty nodal points (joints). Although the mesh is quite 
coarse, it adequately describes the data input for the given 
boundary situation. 
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Figure 28 shows in detail the input and input format 
required for TRIMEG as described in Appendix C, while 
Figures 29 and 30 present the output and output format. 

Each SEL contains one row, three columns and three 
slices . 
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26 



27 



28 
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V 



Figure 22. Crack Propagation Problem. 
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Figure 24. Crack Propagation Problem. 
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Cartesian SEL Diagram 



Figure 25. Folded Plate. 
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Figure 26. Local SEL Diagram. 
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79 




Final Mesh 



Figure 27. Folded Plate. 
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MESH PARAMETERS 
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Figure 30 (continued) 
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700 NSTCP= 1 

V>R I TE (6,725) MAXKEL 
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SUPER ELEMENT COLUMN LOOP 
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25C CC 2b0 N=1 »4 
26C XI(N)=XIVAL 

CC 270 N = 5 » 12 , 7 
2 7 C XI(M=XIVAL+ OELXI 
CC 280 N = 6 t L 1 » 5 
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V 1 = XYC ( I , 2 ) 

500 V 4 L( I) = FCC(X,Y,X 1 ,Y 1 ) 
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XMAX=DMAXl(XMAX»CGRD ( 1 , 1 )) 
YMAX=DKAX 1 (YMAX»CCRD (I , 2 ) ) 
XNIN=OMIN 1 (XMIN,CORD(I, 1 ) ) 



VMN=DMIN 1 (YMIN,CCRD(I , 2 ) ) 
IF (XMIN.GE.O.DO) GO TO 40 
XPAX-XNAX-XMIN 
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APPENDIX F 



COMPUTER LISTING (TRIMEG- 1) 
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IF (LLS .EC .LLS2 ) NN=N+4 
202 l\C(M=NK 
N1=NC(1) 

N 2 = NC ( 2 ) 
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IF ( (MS(7) .EQ.l) .AND. ( MS ( 6 ) .EQ.OJ ) NC1=NC1+1 

IF ( ( (MS(7 ) .EQ.O) . AND. IMS(8) .EQ.l)) .AND. ( P F (K A , 6 ) . EQ . 1 ) ) NC1 = NC1+1 
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1( .SD1*( X*X+Y*Y+Z*Z)-1.9C1)/6.4D1 
FCM(X>Y»Z»X1»Y1»ZL) = ( 1 . OD C-X*X ) * ( 1 . OD 0+ . 9 C 1*X*X 1 ) * ( L . ODO + Y * Y 1 ) * 
1 ( 1.CD0 + Z-Z1 )*.9D1/6.4DI 
IF (NBNODE .EQ .32 ) GC TO 500 
IF { NBNCDE .EQ .20 ) GO TO 200 
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400 V£L(I)=FQMIZ,X,Y,XI,Y1) 
RETURN 
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APPENDIX G 



COMPUTER LISTING (TRIMEG-2) 
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IF (NSTCP.EC.l) GC TO 5 
CALL COORD 

IF (NPLCT.NE.l) GO TO 5 
CALL TRFR(TETA, ALPHA, BETA) 
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IF ( ( MS ( 7 ) .EQ.L) .ANC.lMS(o) .EQ.OJJ NCI=NC1+1 

IF ( ( l MS (7 i .EQ.O ) .AND. (MS (8 ) .EQ.L)) . AND . ( N F ( K A , 6 ) .EQ.L) ) NC1=NC1+I 
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APPENDIX H 



COMPUTER LISTING (TRIMEG- 3) 
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CC 60 N = 1 ,NSEL 

IF (ISiPTSBT .GT.MAXBJT) NST 0 P =2 
R E AC ( 5 , 11 ) I, (NSCONt I , J ) ,J=i, 12 ) 
WRITE ( 6 , 31 ) I , (NSCCM I , J ) ,J = 1 , 13 ) 
NPTS 6 T = NPTSBT + NSC 0 i\( 1 , 13 ) 
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K l=K+l 
I 1 = 1 + 1 
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165 IF ( II .EC. 1) GO TO 300 

IF ( ( MS ( 2 ) .EQ.O) .AND. (MS ( 3 ) . EQ • 0 ) ) GO TO 3CC 
IF (M.NE.ML) GO TO 171 
M V 2 = MEL(KL ,3) 
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GC TO 243 

236 IF < < (M .EQ.M1) .AND . (MS(2) .EQ. 1) ) .AND. (MFE(KL,4) .EQ.l )) GO TO 238 
GC TO 240 

238 NCCN(MBS3» L) = AC0MNBS6,21) 

NCCNC M, 4)=NC0M MBS3 , 1 ) 
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256 IF (M.NE.KL) GO TC 258 

IF <{MF(K, l) .EG.i) .AND. <MF(K»8) .EC.l) ) GC TC 257 
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